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1 Introduction 

To arrive at a Hilbert space representation of measurements starting from plausible "phys- 
ical" axioms has been a dream of many physicists and mathematicians for almost seventy 
years. Of course, one could not expect to recognize the axioms from nothing but exper- 
imental data because they — provided they exist — must be rather involved. Therefore the 
scientists took the opposite road by starting with Hilbert space and trying to read off essen- 
tial mathematical properties so as to be able to eventually simplify them and arrive at the 
simple physically plausible axioms. 

The first breakthrough along the opposite-road was made by Birkhoff and von Neumann 
in 1936 m who recognized that a modular lattice which can be given a physical background 
underlies every finite dimensional Hilbert space. In the early sixties Mackey [^ (and Zierler 
01) arrived at six axioms for a poset (partially ordered set) of physical observables which 
he essentially read off' from the Hilbert space properties. In an additional famous seventh 
axiom, he then postulated that the latter poset be isomorphic to the one of the subspaces 
of an infinite-dimensional Hilbert space. A few years later PironQ], MacLaren ||^, Amemyia 
and Araki 0, and M§czyhski |^ starting with such a poset with infima and suprema of 
every two-element subset (lattice) and using similar axioms, proved that the lattice — usually 
called the Hilbert lattice — is isomorphic to a pre-Hilbert space. That enabled M§,czyhski [0 
to postulate only the kind of a field over which the Hilbert space should be formulated: he 
chose the complex one. The Hilbert lattice is then a lattice of the subspaces of the Hilbert 
space. 

At that time it seemed that only two other fields could have been postulated: the real 
and the quaternionic ones. But in the early eighties Keller p] showed that there are other 
non-standard (non-archimedean) fields over which a Hilbert space can be defined. Also, the 
axioms themselves proved to be too complicated to be given plausible physical support or 
simplified. Thus, the whole project lost its appeal and the majority of researchers left the 
field. However, in 1995 Maria Pia Soler |^ proved that an infinite dimensional Hilbert space 
can only be defined over either real, or complex, or quaternionic field (i.e., that only finite 
dimensional ones allow non-standard fields). 

The latter result renewed the interest in the problem of reconstructing the Hilbert space 
from an algebra of observables. [|10], 0, |T2|, |13|, 0, [l^ Also, recently devised quantum com- 



puters prompt for such a reconstruction from an algebra which is in the field of quantum 
computing usually called quantum logic in analogy to classical logic underlying classical com- 
puters. In particular, if we wanted quantum computers to function as quantum simulators, 
i.e., to directly simulate quantum systems through their description in the Hilbert space, 
we apparently have to start from such an algebra. For, this would be the only presently 
conceivable way of typing in the Hamiltonian at the console of the quantum simulator. This 
also means that we have to go around the present standard axioms for the Hilbert lattice not 
any more because they are complicated and physically non-grounded but because they in- 
clude universal and existential quantifiers which are unmanageable by a quantum computer. 
A way to do so would be to find lattice equations as substitutes for the axioms. Hilbert 
lattice satisfies not only the orthomodularity equation but a number of other equations as 
well. Thus, if we started with a such a lattice — which can easily be physically supported 



by e.g. a quantum computer design — we could obviously simplify the axioms and possibly 
ultimately dispense with them. The problem is that only two groups of the equations satis- 
fied by Hilbert lattices — i.e., in any Hilbert space — have been found so far. We do not know 
whether the Hilbert space equations form a recursively enumerable set, i.e., whether we can 
determine them all. What we can do however is to try to find as many such equations as 
possible and group them according to their recursive algorithms. Each such equation can 
simplify the present axioms. 

However, since already equations with 4 variables contain at least about 30 terms which 
one cannot further simplify, a proper tool for finding and handling the equations is indis- 
pensable. As a great help came Greechie diagrams (condensed Hasse diagrams) which we 
will define precisely later on. E.g., to find that two equations cannot be inferred from each 
other it suffices to find two Greechie lattices which the equations interchangeably pass and 
fail. As an illustration of how "easily" one can find a lattice without a computer program we 
cite Greechie himself: "[In 1969] a student, beginning his dissertation, found such a lattice. 
It was terribly complicated and had about eighty atoms. The student left school and the 



example was lost. I've been looking for one ever since. Recently [in 1977!] I found one." ^Q\ 
So, we need an algorithm for finding Greechie diagrams and another for finding whether a 
particular equation passes or fails them. In this paper we give both. They would not only 
support the afore-mentioned project of obtaining the Hilbert space from physically plausible 
axioms but would also serve for obtaining new equations in the theory of Hilbert spaces. 

The first attempt at automated generation of Greechie diagrams was made in the early 
eighties by G. Beuttenmiiller a former student of G. Kalmbach. |1^, pp. 319-328] The 
algorithm itself is not given in the book but G. Beuttenmiiller kindly sent us the listing of 
its translation into Algol. We rewrote it in C, and with a fast PC it took about 27 days to 
generate Greechie diagrams with 13 blocks. We estimated it would take around a year for 
14 blocks and half a century for 15 blocks, so we looked for another approach. 



The technique of isomorph-free exhaustive generation |I^ of Greechie diagrams gave us 
not only a tremendous speed gain — 48 seconds, 6 minutes, 51 minutes, 8 hours and 122 
hours for 13-17 blocks, respectively (for a PC running at 800 MHz) — but also essentially 
new results: Beuttenmiiller's algorithm must be at least incomplete since the numbers of 
non- isomorphic Greechie diagrams in Kalmbach's book {^^, p. 322) are wrong. In Sec. ^ 
we give the algorithm for the above generation. 

In Sec. ^ we give an algorithm for checking whether a particular equation fails or passes 
in Greechie diagrams provided by the algorithm from Sec. 0. The algorithm has helped us 



to find new equations that hold in any infinite dimensional Hilbert space. [19 



2 Isomorph-free exhaustive generation of Greechie di- 
agrams 

The following definitions and theorem we take over from Kalmbach [O] and Svozil and 



Tkadlec [^. Definitions in the framework of quantum logics (a-orthomodular posets) the 
reader can find in the book of Ptak and Pulmannova. 



Definition 2.1. A diagram is a pair {V, E), where V ^ ^ is a set of atoms (drawn as points) 
and E C exp V \ {0} is a set of blocks (drawn as line segments connecting corresponding 
points). A loop of order n > 2 (n being a natural number) in a diagram {V, E) is a sequence 
(ei,...eb) G -E" of mutually different blocks such that there are mutually distinct atoms 
i/i, . . . , z/„ with UieCiH Ci+i {i = l,...,n, Cn+i = Ci). 

Definition 2.2. A Greecliie diagram is a diagram satisfying the following conditions: 

(1) Every atom belongs to at least one block. 

(2) If there are at least two atoms then every block is at least 2- element. 

(3) Every block which intersects with another block is at least 3-element. 

(4) Every pair of different blocks intersects in at most one atom. 

(5) There is no loop of order 3. 

Theorem 2.3. For every Greechie diagram with only finite blocks there is exactly one (up to 
an isomorphism) orthomodular poset such that there are one-to-one correspondences between 
atoms and atoms and between blocks and blocks which preserve incidence relations. The poset 
is a lattice if and only if the Greechie diagram has no loops of order 4- 

In the literature, a block is also called an edge and an atom is also called a vertex or 
node. (However, we reserve the term node for an element of a Hasse diagram.) 

From the above definitions it is clear that a block can have not only 3 atoms but also 2 or 4 
or more atoms. However, practically all examples of Greechie diagrams used in lattice theory 
are nothing but pasted 3-atom blocks. We are aware of only two important contributions 
containing 4-atom blocks (two proofs of the existence of finite lattices admitting no states 
given in |2l|. Fig. 2.4.5, p. 37] and [|17[ Fig. 17.3, p. 275] and of only one result for n and oo 



giving orthomodular lattices without states |T^, Fig. 17.4, p. 275]. 



For this reason, we have initially focussed on generation of diagrams with every block 
having size 3. Nevertheless, our description of the generation algorithm will allow larger 
blocks in anticipation of the next version of our generation program. Our program for 
checking equations already handles large blocks. Also, we are interested only in the diagrams 
which correspond to lattices, i.e., only in those containing no loops of order 4. Since the 



condition (5) of Def. p.2| states that there are no loops of order 3, this means that we 
are interested only in diagrams with loops of order 5 and higher. Those 3-atom Greechie 
diagrams which correspond to lattices we call Greechie- 3-L diagrams. 

A diagram is connected if, for each pair of atoms u, u', there is a sequence of blocks 
61,62, • • • , 6fc such that u G 61, u' G Ck and Ci fl 6i+i 7^ for 1 < i < k—1. In Section ^ 
we will illustrate how the properties of unconnected diagrams are not necessarily a simple 
combination of the properties of their connected components, so our algorithms will handle 
both connected and unconnected diagrams. An isomorphism from a diagram {Vi,Ei) to a 
diagram (V2, E2) is a bijection from Vi to V2 such that induces a bijection from Ei to E2. 
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The isomorphisms from a diagram D to itself are its automorphisms, and together comprise 
its automorphism group Aut(D). 

\i D = {y, E) is a diagram and e G £", then D — e is the diagram obtained from D by 
removing e and also removing any atoms that were in e but in no other block. Conversely, 
if e is a set of atoms (not necessarily all of them atoms of D), then D + e is the diagram 
{yUe.EU {e}). Clearly {D + e)-e = D. 

We will describe the generation algorithm in some generality to assist future applications. 
Suppose that C is some class of diagrams closed under isomorphisms (for example, connected 
Greechie-3-L diagrams). \i D = (V, E) E C and \E\ > 1, there may be some e E E such that 
D — e E C. If there is no such block, we call D irreducible. It is obvious that all diagrams in 
C can be made from the irreducible diagrams in C by adding a sequence of blocks one at a 
time, all the while staying in C. Such a sequence of diagrams is a construction path for D. 

The basic idea behind our algorithm is to prune the set of construction paths until (up 
to isomorphism) each diagram in C has exactly one construction path. This is achieved by 
two techniques acting in consort. 

The first technique is to avoid equivalent extensions. Suppose D E C and ei, 62 are such 
that D + ei,D + 62 E C. D + ci and D + e2 are called equivalent extensions of D ii \ei\ = \e2\ 
and there is an automorphism of D which maps V Cici onto V (162- It is easy to see that the 
equivalence of ei and 62 implies the isomorphism oi D + ci and Z) + 62, by an isomorphism 
that takes ei onto 62, so we do not lose any isomorphism types of diagram if we make only 
one of them. 

The second technique is somewhat more complicated. Suppose we have a function m 
with the following properties. 

• m{ ) takes a single argument D which is a reducible diagram in C. It returns a value 
which is an orbit of blocks under the action of Aut{D). 

• D — e E C for every e G m{D). 

• If D' is a diagram isomorphic to D, then there is an isomorphism from D to D' that 
maps m{D) onto m{D'). 

We will explain how to compute such a function m{ ) later; for now we will describe its 
purpose. Take any reducible D E C and e G m{D), then form D — e. If we do the same 
starting with a diagram D' isomorphic to D — take e' G m{D'), then form D' — e' — the third 
property of m{ ) implies that D — e and D' — e' are isomorphic. Thus, the function m{ ) 
enables us to define a unique isomorphism class, that oi D — e for e G m{D), as the parent 
class of the isomorphism class of D. Since we wish to avoid making isomorphism types more 
than once, we can decide to only make each diagram from its parent class. If we happen to 
make it from any other class, we will reject it. 



The result of the theory in |]T8| is that the combination of the above two techniques results 
in each isomorphism class being generated exactly once. The precise method of combination 
is given by the following algorithm. 



Definition 2.4. Isomorph-free Greechie diagram generation procedure 

procedure scan [D : diagram; j3 : integer) 

if D has exactly f3 blocks then 
output D 

else 

for each equivalence class of extensions D + e do 
if e G m{D + e) then scan(D + e,(3) 

end procedure 

In Figure [^ we show the top four levels of the generation tree as produced by our 
implementation of the algorithm for connected Greechie-3-L diagrams. The lines joining the 
diagrams show the parent-child relationship. Between a parent and its child, one block is 
added. Note that diagram D43 is made by adding a block to D^i, but it could also be made 
by adding a block to -D3,2- The reason that D31 is its real parent is that mlD^^s) consists of 
the upper right and lower right blocks (which are equivalent) of .04,3. (This is a fact of our 
implementation which cannot be seen by looking at the figure.) When D^^^ is made from 
D^i, the new edge is seen to be in m^D^^^) and so the diagram is accepted. When it is made 
from -03,2; the new edge is found to be not in mlD^^^) and so the diagram is rejected. The 
idea is that each (isomorphism type of) diagram is accepted exactly once, no matter how 
many times it is made. This is proved in the following theorem. 

Theorem 2.5. Suppose we call scan{D, (3) for one D from each isomorphism class of irre- 
ducible diagram in C that has at most j3 blocks. Then the output will consist of one diagram 
from each isomorphism class in C with exactly (3 blocks. 



Proof. The theorem is a special case of one in |]T8[ , and the reader is referred to that paper 



for a strictly formal proof. Here we will give a slightly less formal sketch. 

Let us say that a diagram D is accepted by the algorithm if a call scan{D, (3) occurs. We 
will first prove that at least one member of each isomorphism class of diagram in C with at 
most (3 blocks is accepted. Then we will prove that at most one member of each isomorphism 
class is accepted. These two facts together will obviously imply the truth of the theorem. 

Suppose that the first assertion is false: there is an isomorphism class in C, with at most 
j3 blocks, that is never accepted. Let Z) be a member of such a missing isomorphism class 
which has the least number of blocks. D cannot be irreducible, since all irreducible diagrams 
are accepted explicitly. Thus, we can choose e G m{D) and consider D — e. Since D — e E C 
and D — e has fewer blocks than D, at least one isomorph D' oi D — e is accepted. 

The isomorphism from D — e to D' maps V{D) fie onto some subset of V{D'). Let e' be 
a set of atoms consisting of that subset plus enough new atoms to make e' the same size as e. 
The for loop considers some extension D' + e" equivalent to D' + e', since it considers all 
equivalence classes of extensions. Moreover, since e G m{D) we can infer that e' G m{D' + e') 
and consequently that e" G m{D' + e"). This means that the algorithm will perform the call 
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Figure 1: Generation tree for connected Greechie-3-L diagrams 



scan{D' + e", (3), which is a contradiction as D' + e" is isomorphic to D and the isomorphism 
class of D was supposed to be not accepted at alL This proves that all isomorphism classes 
are accepted at least once. 

Next suppose that some isomorphism type is accepted twice. Namely, there two isomor- 
phic but distinct diagrams D and D' in C, with at most M blocks, such that both D and D' 
are accepted. Choose such a pair D, D' with the least number of blocks. 

As before, D and D' cannot be irreducible, so they must be accepted by some calls 
scan{{D — e) + e, (3) and scan{{D' — e') + e', (3) which arise from the calls scan{D — e, (3) and 
scan{D' — e', (3), respectively, where e G m{D) and e' G m{D'). The properties of r?7,( ) ensure 
that D — e and D' — e' are isomorphic, so they must in fact be the same diagram D" (since 
isomorphism classes with fewer blocks than D are accepted at most once by assumption). 
However, D" + e and D" + e' are equivalent but distinct extensions of D", which violates the 
for loop specification. This contradiction completes the proof. D 

The success of the algorithm requires us to be able to find the irreducible diagrams in C 
by some other method, but in many important cases this is easy. We give the most important 
example. 

Theorem 2.6. Suppose C is a class of Greechie diagrams defined by some fixed set of per- 
missible block sizes, some fixed set of permissible loop lengths, and an optional restriction to 
connected diagrams. Then the only irreducible diagrams in C are those with one block. 

Proof. Consider a diagram D & C with more than one block. 

If C is not restricted to connected diagrams, D — e E C for any e G E{D), so D is 
reducible. 

Suppose instead that C contains only connected diagrams. Choose a longest possible 
sequence 5* of distinct blocks ei, 62, . . . , e^, where e, fl Cj+i 7^ for 1 < z < A;— 1. Let ui and 
1/2 be two atoms oi D — Ck- Since D is connected, there is a chain of blocks from ui to z/2. 
This same chain is in D — e^ unless it contains e^. However, all the blocks D intersecting 
Cfc are in S (or else S can be made longer), so Ck can be replaced in S by some portion of S. 
Hence D — Ck is connected, so D is reducible. D 

The correctness of the algorithm does not depend on the definition of m( ) provided it 
has the properties we required of it. The actual definition of m{ ) used in our program is 
carefully tuned for optimal observed performance, and is too complicated to describe here 
in detail, but we will outline a simpler definition that is the same in essence. 

The key to our implementation of m( ) is the first author's graph isomorphism program 
nauty [^ can be used, nauty takes a simple graph G, perhaps with colored atoms, and 
produces two outputs. One is the automorphism group Avit{G), in the form of a set of 
generators. The other is a canonical labelling of G, which is a graph c{G) isomorphic 
to G. The function c is "canonical" in the sense that c{G) = c{G') for every graph G' 
isomorphic to G. To apply nauty to a diagram {V, E) , we can use the incidence graph 
G = {VUE,{{v,e)\vee}). 



The generators for Aut(G') can be easily converted into generators for Aut(D) and then 
used to determine the equivalence classes of extensions. This enables us to implement the 
requirement of avoiding equivalent extensions. 

The canonical labelling c( ) produced by nauty enables us to define m{ ). Take the block 
e such that D — e & C and e is given the least new label by c( ). If we define m{D) to be the 
orbit of blocks that contains e, we find that the three requirements we imposed on m{ ) are 
satisfied. 

As we have said, our real program uses a more complex definition of m{ ). We do not 
use the incidence graph G, but instead use a prototype variant of nauty that operates on 
diagrams directly. Since our program makes connected diagrams, we took m{D) to be an 
orbit of feet if there were any, where a foot is a block with only one atom that also lies 
in other blocks. This avoids many connectivity tests, since removal of a foot necessarily 
preserves connectivity. It also avoids many futile extensions: adding a non-foot e must be 
done in such a way that any existing feet become non-feet, as otherwise e ^ 171(0 + e). 

In order to generate only the connected Greechie-3-L diagrams having M blocks but no 
feet, a reasonable approach is to generate all the diagrams, having feet or not, with M — 1 
blocks first. Then the M-th block can be added in such a way that uses at least 2 of the 
existing atoms and also turns any feet into non-feet. It is also possible to make a generator 
that makes foot-free diagrams while staying entirely within that class, but it does not appear 
likely to be much different in efficiency. 

Program greechie. Our implementation of the algorithm is a self-contained program called 
greechieQ that takes as parameters the number of blocks, an optional upper bound on the 
number of atoms, and whether or not feet are permitted. It then produces one representative 
of each isomorphism class of connected Greechie-3-L diagram with those properties. The 
diagrams can then be processed as they are generated, with no need to store them. There is 
also an option for dividing the set of diagrams into disjoint subsets, and efficiently producing 
only one of the subsets. This allows long computations to be broken into manageable pieces 
that can be run independently, even on different computers, without much change to the 
total running time. 

In Table |l| we list the numbers of Greechie-3-L diagrams for small values of a and f3. 
In each cell of the tables, the upper value is the total number of connected Greechie-3-L 
diagrams, and the lower value is the number of those which have no feet. Both counts are 
if the table cell is empty. The table includes all possible values of /3 for a < 29 and all 
possible values of a for /3 < 17. 

For reasons explained later, we have particular interest in those diagrams containing 
close to the maximum number of blocks for a given number of atoms. This prompted us to 
compute additional near-maximal diagrams past the size where finding all the diagrams is 
practical. 

To keep the discussion simple, we restrict ourselves to Greechie-3-L diagrams, not neces- 
sarily connected. By the type of a diagram we mean the pair {a,f3), where a is the number 



[ftp : //m3k . grad . hr/pavicic/greechie , attp: //cs . anu. edu. au7|~bdm/nauty/greechie .html. Many 



of the diagrams computed with the program are also available at those places. 
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1448834695 
4690081789 


6661929716 
10298834720 


7447274324 
8422315646 


1916771296 
1956383755 


82563871 
82670819 


359550 
359550 


245 
245 


17618049369 
26495317590 



Table 1: Counts of connected Greechie-3-L diagrams 
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of atoms and (3 is the number of blocks. The rank of an atom is the number of blocks which 
contain it. 

The first observation is that a diagram D of type (a, f3) has an atom whose rank is at 
most [3/5/aJ . Over most of our computational range, a < /?, so this value is at most 2. If 
there is an atom of rank 1, we can make D either by adding a foot to a diagram of type 
(a— 2,/3— 1), or by adding one block and one atom to a diagram of type (a— I,/?— 1). On 
the other hand, if there is an atom of rank 2 but none of rank 1, we can make D by adding 
one atom and two blocks to a diagram of type (a— 1, (3—2). So, if we have already made the 
diagrams of types (a— 2,/5— 1), (a— 1,/3— 1), and (a— 1,/5— 2), we can easily extend them to 
make those of type (a, (3), provided a < 13. 

Sometimes the class (a— 1,/5— 2) may be too onerous to compute. In this case, there 
are two other approaches we might be able to take to making those diagrams of type (a, /3) 
whose minimum atom rank is 2. Define a^ to be the number of atoms of rank i (and recall 
that we are assuming ai = 0). Counting the pairs (block, atom in block) in two ways, we 
have 2a2 + 3^3 + 4(a; — 0:2 — a^) < 3/3. Since also 0^2 + «3 < a, we have 



2^2 + as > 4a — 3/3 and a2 > 3a — 3/3. 



fll 



Now consider the case of a Greechie-3-L diagram with ai = and 6a > 7/3. Applying 
the second part of (|l]) we find that 2a2 > /3, which implies that some block contains at least 
two atoms of rank 2. Therefore, we can make the diagram by adding two atoms and three 
blocks to a diagram of type (a— 2,/?— 3). 

If 6a < 7/3 but a > /3, there might be no atoms of rank 1, nor two atoms of rank 2 in 
the same block. In this case, we know that there are exactly 2a2 blocks containing an atom 
of rank 2. The total rank of the atoms of rank greater than 3 is 3/3 — 2a2 — 3a3, so we have 
at least 4a2 — (3/3 — 2a2 — 3a3) = 6a2 + 3a3 — 3/3 pairs (x, y) such that x and y are atoms of 
rank 2 and 3, respectively, lying in the same block. Now, if we suppose that 14a > 15/3, we 
find that 6a2 + 3a3 — 3/3 > a^, implying that two of the pairs (x, y) have the same y. That 
is, there is an atom of rank 3 lying in two blocks which each contain an atom of rank 2. 
Therefore, we can make this diagram by adding 3 atoms and 5 blocks to a diagram of type 
(a-3,/3-5). 



a\(3 


25 


26 


27 


28 


29 


30 


31 


32 


33 34 


35 36 


30 


3982 


4 


















31 


? 


81068 


71 


1 














32 


? 


? 


>313813 


1643 














33 


? 


? 


? 


? 


>51643 


66 










34 


? 


? 


? 


? 


? 


>185733 


2113 


19 






35 


? 


? 


? 


? 


? 


? 


? 


>70035 


325 17 


5 


36 


? 


? 


? 


? 


? 


? 


? 


? 


? >7871 


136 1 


37 


? 


? 


? 


? 


? 


? 


? 


? 


? ? 


? >1693 



Table 2: Counts of large extremal Greechie-3-L diagrams 
Using these ideas, we were able to compute (see Table ||) the Greechie-3-L diagrams with 
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up to 36 atoms having the maximuin and, in some cases, near-maximum possible numbers 
of blocks. 

The maximum number of blocks shown in each row of the table is the maximum possible. 
All the diagrams counted in the table turned out to be connected, even though our programs 
did not assume connectivity. 




Figure 2: The unique Greechie-3-L diagram of type (36, 36) 



Of some interest is that the Greechie-3-L diagrams of type (35, 35) and (36, 36) contain 
only atoms of rank 3. There is a connection here to known results in graph theory, as follows. 
Let D he a diagram in one of these two classes. Define a graph G whose vertices are the 
atoms and blocks oi D. A vertex which is an atom is adjacent in G to a vertex which is 
a block provided the atom lies in the block, and there are no other edges. (This is the 
same incidence graph we defined earlier.) Then each vertex of the graph has valence 3, and 
moreover there are no cycles of length 9 or less (which is equivalent to the requirement that 
D has no loops of length 4 or less). In graph theoretic language, G is a bipartite cubic (or 
trivalent) graph of girth at least 10. The smallest such graphs are three with 70 vertices 
p3| . In two of the three graphs, the functions of "atom" and "block" can be interchanged to 
produce different Greechie diagrams, but in the third this interchange gives an isomorphic 
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diagram. That is how the three graphs correspond to the five diagrams we found. 

Figure shows the unique Greechie-3-L diagram of type (36, 36). Using a shght adapta- 



tion of the program described in [^], we have found that there are no Greechie-3-L diagrams 
of type (37, 37) with every atom having rank 3, and exactly eight such diagrams of type 
(38,38). 



3 Testing conjectures with Greechie diagrams 

In order to test equations conjectured to hold in various classes of orthomodular lattices, it 
is useful to automate their checking against Greechie diagrams. This will let us either falsify 
the conjecture or give us some confidence that might hold in the class of interest before we 
attempt to prove it. Also, finding a lattice in which one equation holds but a second one 
fails gives us a proof that the second is independent. 

To this end we use a program called latticegQ which will check to see if an equation or 
inference holds in each of the Greechie diagrams in a list provided by the user. 

The latticeg program internally converts a Greechie diagram to its corresponding Hasse 
diagram and tests all possible assignments of nodes in the Hasse diagram to the equation. 
In general Greechie diagrams correspond to Boolean algebras "pasted" together. 

1 




x' X 



X 

o 



X 





Figure 3: Greechie diagrams for Boolean lattices 2^, 2'^, and 2^, labeled with the atoms of 
their corresponding Hasse diagrams shown above them. (2^ was adapted from [^, Fig. 18, 
p. 84].) 

The Hasse diagrams for the Boolean algebras corresponding to 2-, 3-, and 4-atom blocks 
are shown in Fig. 0. The Greechie diagram for a given lattice may be drawn in several 



'"^Available at ftp : //ftp . shore . net/members/ndm/quantum-logic as the ANSI C program latticeg. c. 
The program is simple to use and self-explanatory with the — help option. A related program lattice. c 
handles general Hasse diagrams and has built-in those lattices we have found most useful for preliminary 
testing of conjectures. These two programs, along with beran.c for computing the canonical orthomod- 
ular form of any two-variable expression, are the primary computational tools we have used for studying 
orthomodular and Hilbcrt lattice equations. 
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equivalent ways: Fig. ^ shows the same Greechie diagram drawn in two different ways, along 
with the corresponding Hasse diagram. From the definitions we see that the ordering of the 
atoms on a block does not matter, and we may also draw blocks using arcs as well as straight 
lines as long as the blocks remain clearly distinguishable. 

1 



X 




w 




Figure 4: Two different ways of drawing the same Greechie diagram, and its corresponding 
Hasse diagram. 

Recall that a poset (partially ordered set) is a set with an associated ordering relation that 
is reflexive {a < a), antisymmetric [a < b,b < a imply a = b), and transitive {a < b,b < c 
imply a < c). An orthoposet is a poset with lower and upper bounds and 1 and an 
operation ' satisfying (i) if a < 6 then b' < a'; (ii) a" = a; and (iii) the infimum ad a' and 
the supremum aU a' exist and are and 1 respectively. A lattice is a poset in which any 
two elements have an infimum and a supremum. An orthoposet is orthomodular ii a < b 
implies (i) the supremum a U 6' exists and (ii) a U (a' fl 6) = 6. A lattice is orthomodular 
if it is also an orthomodular poset. For example. Boolean algebras such as those of Fig. ^ 
are orthomodular lattices. A a -orthomodular poset is an orthomodular poset in which every 
countable subset of elements has a supremum. An atom of an orthoposet is an element a 7^ 
such that b < a implies 6 = 0. 

In the literature, there are several different definitions of a Greechie diagram. For exam- 
ple, Beran ([^ p. 144]) forbids 2-atom blocks. Kalmbach (|jl^, p. 42]) as well as Ptak and 
Pulmannova |21|, p. 32] include all diagrams with 2-atom blocks connected to other blocks 
as long as the resulting pasting corresponds to an orthoposet. However, the case of 2-atom 
blocks connected to other blocks is somewhat complicated; for example, the definition of a 



loop in Def. |2T^ must be modified (e.g. [11^, p. 42]) and no longer corresponds to the simple 



geometry of a drawing of the diagram. The definition of a Greechie diagram also becomes 
more complicated; for example a pentagon (or any n-gon with an odd number of sides) made 
out of 2-atom blocks is not a Greechie diagram (i.e. does not correspond to any orthoposet). 

The definition of Svozil and Tkadlec [^ that we adopt, Def. |2]^, excludes 2-atom blocks 
connected to other blocks. It turns out that all orthomodular posets representable by Kalm- 
bach's definition can be represented with the diagrams allowed by Svozil and Tkadlec's 
definition. But the latter definition eliminates the special treatment of 2-atom blocks con- 
nected to other blocks and in particular simplifies any computer program designed to process 
Greechie diagrams. 

Svozil and Tkadlec's definition further restricts Greechie diagrams to those diagrams rep- 
resenting orthoposets that are orthomodular by forbidding loops of order less than 4, unlike 
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the definitions of Beran and Kalmbacfi. Tfie advantage appears to be mainly for convenience, 
as we obtain only those Greechie diagrams that correspond to what are sometimes called 
"quantum logics" (a-orthomodular posets). (We note that the term "quantum logic" is 
also used to denote a propositional calculus based on orthomodular or weakly orthomodular 
lattices. 0) 

The definition allows for Greechie diagrams whose blocks are not connected. In Fig. |^ we 
show the Greechie diagram for the Chinese lantern M02 using unconnected 2-atom blocks. 
This example also illustrates that even when the blocks are unconnected the properties of the 
resulting orthoposet are not just a simple combination of the properties of their components 
(as one might naively suppose), because we are adding disjoint sets of incomparable nodes to 
the orthoposet. As is well-known ([l^, p. 16]), M02 is not distributive, unlike the Boolean 
blocks it is built from. 



x' X y y 

O^ O^ 

(a) 


Figure 5: Greechie diagram for the lattice M02 and its Hasse diagram. The dashed line 

indicates that the unconnected blocks belong to the same Greechie diagram. 

The latticeg program takes, as its inputs, a Greechie diagram ASCII representation 
(or more precisely a collection of them) and an equation (or inference) to be tested. This 
ASCII representation is compatible with the output of the programs described in Section |[ 
Currently latticeg is designed to work only with Greechie diagrams corresponding to lat- 
tices, i.e. that have no loops of order 4 or less, as these are the most interesting for studying 
equations valid in all Hilbert lattices. It converts the Greechie diagram to its corresponding 
Hasse diagram (internally stored as truth tables). Finally, the program tests all possible 
assignments of nodes in the Hasse diagram to the input equation under test. 

The latticeg program also incorporates classical propositional metalogic and predicate 
calculus to allow the study of such characteristics as atomicity and superposition. 

As a simple example of the operation of the latticeg program, we show how it verifies 
the passage and failure of the modular law [Equation ([Tl|) below] on the lattices of Figs. |^ 
and 0c. We create a file with a name such as m.gre to represent the lattices, containing the 
lines 

123,345. 
123,345,567. 

and run the program by typing 

latticeg -i m.gre " (av(b"(avc)))=((avb)" (avc)) " 

The program responds with 
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The input file has 2 lattices. 

Passed #1 (5/2/12) 

FAILED #2 (7/3/16) at (av(f ~ (avb) ))=( (avf ) " (avb) ) 

The notation should be more or less apparent but is described in detail by the program's 
help. The numbers "5/2/12" show the atom/block/node count, and the failure shows the 
internal Hasse diagram's nodal assignment to the equation's variables. 

Let us consider an application of the program. Closed subspaces Tia^'Hh of any infinite 
dimensional Hilbert space Ti form a lattice in which the operations are defined in the following 
way: a' = Tif, aHb = TCaClTi-b, and aUb = {Ha + Tih)'^^- In such a lattice its 
elements satisfy the following condition, i.e., in any infinite dimensional Hilbert space its 
closed subspaces satisfy the following equation which is called the orthoarguesian equation: 

a ±b & c ± d & e ± f =^ 

(a U 6) n (c U rf) n (e U /) < 
6 U (a n (c U (((a U c) n (6 U d)) n (((a U e) n (6 U /)) U ((c U e) H (rf U /)))))) (2) 

where a _L 6 means a <b' . 

We wanted first, to reduce the number of variables in this equation and second, to 
generalize the equation to n variables. In attacking the first problem the program helped 
us to quickly eliminate dead ends: a failure of a conjectured equation in a lattice in which 
Equation |^ held meant that the latter equation was weaker and vice- versa. Thus we arrived 
at the following 4-variable equation — which we call the 40A law: 

(4) 

(ai ^1 as) n (01 = 02) < 02 ^1 as • (3) 

(4) 

where the operation = is defined as follows: 



where 



(^) def , (3) , ,, (3) ,^, (3) 

ai = a2 = (01 = 02) U ((01 = 04) n (02 = 04)) , (4) 



oi=02 = ((oi ^1 03) n (02 ->i 03)) u ((o'l -^1 03) n (02 ^1 03)), (5) 



where o — >i b = o' U (o fl 6). We then proved "by hand" that Equations (Ep and (0) are 



equivalent. |]T9l We also proved that the following generalization (which we call the nOA 
law) of Equation (||) 



where 



(oi ^1 03) n (oi =02) < 02 ^1 03 (6) 



Oi = 02 = (oi = 02) U ((oi = o„) n (02 = an)), n>4 (7) 



holds in any Hilbert lattice. To show that this generalization is a nontrivial one, the program 
is all we need — we need not prove anything "by hand." It suffices to find a Greechie diagram 
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Figure 6: Greechie diagram for OML L46. 



in which the 40A law holds and 50A law fails. An 800 MHz PC took a few days to find 
such a lattice (shown in Fig. ^, [r^| while to find it "by hand" is, due to the number of 
variables in the equation and nodes in the corresponding Hasse diagram, apparently humanly 
impossible. 

Considerable effort was put into making the program run fast, with methods such as 
exiting an evaluation early when a hypothesis of an inference fails. Truth tables for all 
built-in compound operations (such as the various quantum implications and the quantum 
biconditional) are precomputed. The innermost loop (which evaluates an assignment) was 
optimized for the fastest runtime we could achieve. The number of assignments of lattice 
nodes to equation variables that must be tested is rt" where n is the number of nodes in 
the Hasse diagram and v is the number of variables in the equation. The algorithm requires 
a time approximately proportional to kii" where k is the length of the equation expressed 
in Polish notation. For a typical equation {k = 20) with no hypotheses, the algorithm 
currently evaluates around a million assignments per second on an 800-MHz PC. The speed 
is typically faster when hypotheses are present, particularly if they have fewer variables than 
the conclusion. If a lattice violates an equation, it often happens (with luck) that the first 
failure will be found quickly, in which case further evaluations do not have to be done. 

The propositional metalogic feature of latticeg, when carefully used in a series of 
hypotheses with a successively increasing number of variables, can sometimes be exploited 
to achieve orders of magnitude speed-up with certain equations containing many variables. 
For example, we tested the 8-variable Godowski equation against a 42-node lattice in 16 
hours, whereas without the speed-up it would have required around 10^ hours. To illustrate 
how this speed-up works, we can add to the 4- variable Godowski equation [Equation ^ below] 
hypotheses as follows: 



{d — *•! a < a — >i d) 



h ~ ((c — >i d) n {d ^1 a) < a — *i d) = 

(a -^1 b) n (6 ^1 c) n (c -^1 d) n (d ^i a) < a 



>i d . 



Here ~ means metalogical not. The hypotheses are redundant in any ortholattice, which 
is the case for the Greechie diagrams that are of interest to us. We take advantage of the 
latticeg feature that exits the evaluation of a lattice nodal assignment if a hypothesis fails. 
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If an assignment to the first hypothesis, with only 2 variables, fails (as it typically does for 
some assignments) it means we don't have to scan the remaining variables. Similarly, if the 
first hypothesis passes but the second (with 3 variables) fails, we can skip the evaluation of 
the 4- variable conclusion. 

We are currently exploring the exploitation of possible symmetries inherent in the Gree- 
chie diagram to speed up the program further, but we have not yet achieved any results in 
this direction. For certain special cases such as the Godowski equations we are also exploring 
the use of a "dynamic programming" technique that may provide a run time proportional 
to kn^ instead of kn^, regardless of the number of variables. 

For fastest run time, it is desirable to screen equations with the smallest Greechie di- 
agrams first. For this purpose what matters is the size of the Hasse diagram and not the 
Greechie diagram. In a chain of blocks each having two atoms connected, a 3-atom block 
adds 4 nodes to the Hasse diagram whereas a 4-atom block adds 12 nodes. For example, 
the decagon (10 blocks) has 42 nodes with 3-atom blocks and 122 nodes with 4-atom blocks. 
With 3-atom blocks, a 6-variable equation — our practical upper limit when there are no 
strong hypotheses — must be evaluated 42^ = 5.5 billion times (a few hours of CPU time 
on an 800-MHz PC), but with 4-atom blocks it would take 122^ = 3.3 trillion evaluations, 
which is currently impractical. So far most of our work has been done using diagrams with 
every block having size 3. 

Two other heuristics have helped us to falsify conjectures more quickly. The first is to 
first scan Greechie diagrams with the highest block-to- atom ratio (Table [^). Such diagrams 
seem to have the most complex "structure" with the most likelihood of violating a non- 
orthomodular equation ( "non-orthomodular equation" here means an equation which turns 
the orthomodular lattice variety into a smaller one when added to it). A drawback is that 
at higher atom counts, virtually every such diagram violates almost any non-orthomodular 
equation, making it very useful for identifying non-orthomodular properties but less useful 
for proving independence results. For example, for 35 and 36 atoms (the case we elaborated 
in Section |^) the highest ratio is 1 (Table ^ and in those diagrams all equations that we 
know to be non-orthomodular fail. Hence, for example 36 x 36 (36 atoms, 36 blocks) is a very 
useful tool for an initial scanning of equations we want to check for "non-orthomodularity." 

The second heuristic is our empirical observation that Greechie diagrams without feet 
often behave the same as the same diagram with feet added. For example, the Peterson 
OML (Fig. 0a), with 32 nodes, is the smallest lattice that violates Godowski's 4- variable 
strong state law [|l9i which holds in any infinite dimensional Hilbert space: 



(a -^1 b) n (6 -^1 c) n (c ^1 d) n {d -^i a) < a ^i d (9) 

but not Godowski's 3-variable law (which also holds in infinite dimensional Hilbert spaces) 

(a -^1 h) n {h —>i c) n (c ^1 a) < a—>ic. (10) 

The Peterson OML is useful as a test for an equation derived from the 4-variable law and 
conjectured to be equivalent to it: if it does not violate the conjectured equation we know 
the equation is weaker the 4-variable law. Now, we observe empirically that we may add a 
foot (a 3-atom block connected at only one point) to any of its 15 atoms without changing 
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this behavior. We have also not seen a chain of feet or combination of feet that changes this 
behavior when added to the diagram. 

So, by scanning only lattices without feet, we can obtain a speedup of 20 times for 14- 
block lattices (Table |^). Supporting this heuristic is the fact that complex Greechie diagrams 
with feet are rarely found in the literature. We obtained an additional support by scanning 
the 40A law given by Equation (^ and Godowski's 3 variable equation (^) through several 
million lattices with free feet, vs. those with feet stripped: we did not find a single difference. 
To our knowledge there is only one special case for which feet do make a difference. Fig. |^ 
shows a lattice that obeys the modular law 



aU (6n (aUc)) 



(a U 6) n (a U c) 



fill 



but violates it when a foot is added (Fig. |^). This special case might well be insignificant 
because of all Greechie diagrams only star-like ones (Fig. ^, D^2 and D^^^^ from Fig. |l|, etc.) 
are modular. As soon as we add any block to any other atom apart from the central one in 
such a lattice we make it non-modular. 





Figure 7: (a) Peterson OML. (b) Greechie diagram obeying modular law. 
diagram violating modular law. 



(c) Greechie 



Another heuristic for reducing the number of diagrams to be scanned is suggested by the 
following observation, although we have not implemented it. We can list the diagrams in such 
a way that (except for the first diagram) each is formed by adding one block to a diagram 
earlier in the list. Whenever the earlier diagram violates an equation, we have observed that 
it is very likely that the new diagram will also violate the equation. By skipping the new 
diagram in this case (presuming the probable violation), a speedup can be obtained. 

However the above observation does not universally hold, i.e. sometimes the new diagram 
will pass an equation violated by the earlier one. An example is shown in Fig. ^j. The diagram 
L38 violates the orthoarguesian law (Equation H). But if we extend L38 by adding two blocks 
as shown in Fig. |p3, the resulting diagram will pass not only this law [equivalent to 40A 
given by Eq. (^] but also 50A and 60A [given by Eq. (H) for n = 5 and n = 6, respectively]. 

If our various speedup heuristics are used for practical reasons, the user must be aware 
that a diagram scan may be incomplete. So if, in our example, the extended L38 (by passing) 
could serve to prove a certain independence result, it would be missed by the scan. Currently 
we have no data to indicate how often such cases would be missed. 
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Figure 8: (a) Greechie diagram for L38; (b) L38 with two blocks added. 

Any scan of diagrams to test an equation is, of course, a priori incomplete since the 
number of diagrams is infinite. The various heuristics we have described may cause some 
lattices in any finite list to be skipped. But if a lattice with the desired properties is found 
more quickly our goal is achieved. Of course a scan can be continued for any diagrams 
omitted by the heuristics if more completeness is desired. 



4 Conclusions 

Greechie diagrams (we used the version given by Definition |2.1| and discussed the others in 
Section ^ are generators of examples and counterexamples of orthomodular non-modular 
lattices. They are special cases of lattices one obtains by using the recent generalization by 
Navara and Rogalewicz p^ E8[ of Dichtl's pasting construction for orthomodular posets and 



lattices |^. Navara and Rogalewicz's method exhaustively generates all finite orthomodular 
non-modular lattices but Greechie diagrams are apparently easier to generate and certainly 
much easier to test lattice equations. For these reasons, Greechie diagrams have been used 
almost exclusively so far. 

Since any infinite dimensional complex Hilbert space is orthoisomorphic to a Hilbert lat- 
tice which is orthomodular and non-modular, Greechie diagrams represent an indispensable 
tool for Hilbert space investigation. This has also been prompted by recent developments in 
the field of quantum computing. However, as we stressed in the Introduction, the existing 
(both manual and automated) constructions of Greechie diagrams and their application to 
Hilbert space properties (which resulted in many important results at the time) recently 
reached the frontiers of human manageability. Therefore, in Section |^ we gave an algorithm 
and a program for generating Greechie lattices with theoretically unlimited numbers of atoms 
and blocks. The algorithm is of a completely different kind from the only earlier algorithm 
and the program is at least 10^ times faster. In Section |^ we then gave an algorithm and 
programs for automated checking of lattice properties on Greechie diagrams. 

Our algorithm for generating Greechie diagrams — given by Definition ^]4| — works by 
defining a unique construction path for each isomorphism class of diagrams. It enabled us 
to produce a self-contained program greechie for automated generation of diagrams with 
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a specified number or range of atoms and/or blocks. Several properties and several special 
types of Greechie diagram construction are discussed in the Section, as are connections to 
some equivalent results in graph theory. 

The algorithm for automated checking of lattice properties described in Section ^ works 
by converting a Greechie diagram to its corresponding Hasse diagram and then converting 
the Hasse diagram to a truth table for supremum and orthocomplementation. It enabled us 
to construct a self-contained program latticeg which takes, as its inputs, Greechie diagrams 
in ASCII representation and an equation or inference or quantified expression to be tested. 
Many programming speed-ups have been used to make the program run as fast as possible. 
For example, all built-in compound operations are precomputed, many C code tricks are 
used, etc. Also, to additionally speed up scanning, we use many Greechie diagram heuristics 
that we found recently: a lattice equation is most likely to fail in a lattice with the highest 
block-to-atom ratio, scanning of an equation on a diagram with no feet and on the same 
diagram with feet added makes no difference for most diagrams, etc. Several properties 
and several special cases of Greechie lattices are given and discussed in the Section. In 
particular, it is explained how we proved that our recent n-variable generalization of the 
orthoarguesian equation is a non-trivial one by using nothing but latticeg applied to an 
output of greechie. 
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